Integrable Flows in 
c = 1 String Theory 



Ivan K. Rostov 

Max- Planck- Institut fur Gravitationsphysik 
Albert- Einstein- Institut 
Am Muehlenberg 1, D-144%4 Golm, Germany 

and 

Service de Physique Theorique, CEA/DSM/SPhT, 
Unite de recherche associee au CNRS, CEA/Saclay 
91191 Gif-sur-Yvette cedex, France] 



In these notes we review the method to construct integrable deformations of the 
compactified c = 1 bosonic string theory by primary fields (momentum or winding modes), 
developed recently in collaboration with S. Alexandrov and V. Kazakov. The method is 
based on the formulation of the string theory as a matrix model. The flows generated by 
either momentum or winding modes (but not both) are integrable and satisfy the Toda 
lattice hierarchy*. 



August, 2002 

° kostov@spht.saclay.cea.fr 

f Permanent address 

* Based on the talks presented at the James H. Simons Workshop on Random Matrix Theory 
Stony Brook, February 20 - 23, 2002, and the International Conference on Theoretical Physics 
Paris, UNESCO, July 22-27, 2002 



Contents 

|L~| Introduction |l| 

|2] The matrix model for the c = 1 string theory revisited |3| 

\2.1. Partition function of the Euclidean theory compactified on a circle with radius R . . |3| 



2.2. | Chiral quantization |4j 

2.3. Gauge-invariant collective excitations: momentum and winding modes |3] 

2.4. The partition function as a three-matrix integral [5] 

|2.5. Integration with respect to the angles |6| 

|3] The partition function of the non-perturbed theory |6| 

|3.1. Eigenfunctions of the non-perturbed Hamiltonian |6| 



3.2. Cut-off prescription, density of states and free energy 



|3.3. Calculation of the free energy using the E'-representation of the operators x± 



fi~| Integrable perturbations by momentum modes 10 

[4.1. | One-particle eigenfunctions and the density of states of deformed system 10 

[4.2.| Dressing operators 11 



4.3. Lax operators and string equation 12 



4.4. Integrable flows 13 



4.5. Representation in terms of a bosonic field 14 



4.6. The dispersionless (quasiclassical) limit 14 



4.7.| One-point correlators in the dispersionless limit |16 



8] Example: sine-Gordon field coupled to 2D gravity 16 



fy\ Integrable perturbations by winding modes 17 



|5.1. The partition function as integral over the U (N) group 17 



5.2. Evaluation of the partition function of the non-perturbed theory 17 



5.3. Toda integrable structure 18 



5.4. Lax operators and string equation 19 



|5.5. Quasiclassical limit and relation to the conformal map problem 19 



|6] Conclusion 20 



|Appendix A. Momentum and winding modes in the c = 1 Euclidean string theory .... 20 



1. Introduction 

The c = 1 string theory (the theory of random surfaces embedded in 1 dimension) 
can be constructed as the collective theory for a 1-dimensional N x N hermitian matrix 
field 0. The U(N) -invariant sector of the matrix model is described by a non-relativistic 
quantum mechanics of free fermions in an upside-down quadratic potential. The mapping 
to free fermions allows to calculate virtually any quantity in the string theory in all orders 
in the genus expansion. 

The elementary excitations of the c = 1 string represent collective excitations of free 
fermions. The tree-level S'-matrix can be extracted by considering the propagation of 
"pulses" along the Fermi sea and their reflection off the "Liouville wall" @. 

The exact nonperturbative S'-matrix has been calculated by Moore et al ||. Each 
S-matrix element can be associated with a single fermionic loop with a number of external 
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lines. One can then expect that the theory is also solvable in a nontrivial, time-dependent 
background generated by a finite tachyonic source. Dijkgraaf, Moore and Plesser |4|] demon- 
strated that this is indeed the case when the allowed momenta form a lattice as in the case 
of the compactified Euclidean theory. In [Q it has been shown that the string theory com- 
pactified at any radius R possesses the integrable structure of the Toda lattice hierarchy 
|J. The operators associated with the momentum modes in the string theory have been 
interpreted in J| as Toda flows. 

The explicit construction of these Toda flows is an interesting and potentially impor- 
tant task, because this would allow us to explore the time-dependent string backgrounds. 
This problem was recently solved in f(|. The method used in || is conceptually similar to 
the method of orthogonal polynomials of Dyson-Mehta in the interpretation of M. Douglas 
0, which has been used to solve the c < 1 matrix models H|J. The construction of 
allows to evaluate the partition function of the string theory in presence of finite perturba- 
tion by momentum modes. The simplest case of such a perturbation is the Sine-Liouville 
string theoryEH. 

Besides the momentum modes, the compactified string theory contains a second type 
of excitations associated with nontrivial windings around the target circle. In the world- 
sheet description, the string theory represents a compactified gaussian field coupled to 2d 
quantum gravity. Then the momentum and winding modes are the electric and magnetic 



operators for this gaussian field. It is known [p!T||12| that the winding modes propagate 
in the non-singlet sector of the c = 1 matrix model, which is no more equivalent to a 
system of fermions. Explicit realization of the winding modes can be given by gauging the 
compactified matrix model. Then the winding modes are realized as the Polyakov loops 



winding around the Euclidean time interval [T^] 



The electric and magnetic operators are exchanged by the usual electric-magnetic 
duality R — > 1/R called also T-duality. If the momentum modes of the compactified string 
theory are described by the Toda integrable flows, then by T-duality the same is true also 
for the winding modes. A direct proof of is presented in ||13|1 , where it is also shown that 
the the grand-canonical partition function of the matrix model perturbed by only winding 
modes is a r-function of the Toda lattice hierarchy. Applying the T-duality backwards, we 
conclude that the same is true for the partition function of the theory perturbed by only 
momentum modes. 

A special case represents the theory is compactified at the self-dual radius R = 1, it 
is equivalent to a topological theory that computes the Euler characteristic of the moduli 
space of Riemann surfaces |14|| . When R = 1 and only in this case, the partition function 



of the string theory has alternative realization as a Kontsevich-type model 0,|T5 



In this notes we summarize the results concerning the Toda integrable structure of 
the compactified c = 1 string theory obtained in refs. []T3| , [T6| , |T7| ,|6|| . We will follow mostly 



the last work |J, but will present the construction in the framework of the Euclidean 
compactified theory, while H discussed the theory in Minkowski space. 



1 It is conjectured [1C] that such a perturbation can lead to a target space geometry with a 
horizon. 



2 



2. The matrix model for the c = 1 string theory revisited 

2.1. Partition function of the Euclidean theory compactified on a circle with radius R 

The c = 1 string theory (see Appendix A) describes the critical behavior of the large- 
N matrix quantum mechanics (see the review |l| and the references therein). The critical 
point is associated with the maximum of the matrix potential. In the scaling limit, which 
is dominated by dense planar graphs, the potential can be approximated by a quadratic 
one (with the wrong sign), which is stabilized by imposing a cut-off wall far from the top. 
The relevant piece of the matrix Hamiltonian is thus 

H Q = ^Tr (P 2 -M 2 ), (2.1) 

where P = —id/dM and My is an iV x N hermitian matrix variable. The cosmological 
constant (jl is introduced as a "chemical potential" coupled to the size of the matrix N, 
which should be considered as a dynamical variable. 

We will consider the Euclidean theory with periodic Euclidean time x = —it. The 
time interval 

P = 2nR (2.2) 

can be also interpreted as inverse temperature. The matrix model describes the c = 1 
string only when (3 > (3kt = 47r. At (3 = (3kt the winding modes associated with 
the iV(iV — 1) angular degrees of freedom becomes important and produce a Berezinski- 
Kosterlitz-Thouless-like phase transition to a c = string theory . The angular degrees 
of freedom can be interpreted as the winding modes of the string, which represent strings 
winding several times around the target circle . 

A sensible matrix model for (3 < (3kt can be constructed by introducing an additional 



SU(N) gauge field A(x), which projects to the sector free of winding modes ||12|1 . The 
partition function of the compactified system is 

ZQm) = f PN ZN> (2-3) 

N>0 

where Zn is given by a functional integral with respect to the one- dimensional matrix 
hermitian fields M{x) and A{x) with periodic boundary conditions 

Z N = J VPVMVA e~^ So . (2.4) 

periodic 

The action functional is 

r/3 



So = Tr / (iPV A M - H ) dx, (2.5) 
Jo 

VaM denotes the covariant time derivative 

V aM = d x M - i[A, M]. (2.6) 

The gauge field can be also used to study perturbations by winding modes, with appro- 
priately tuned coupling constants, so that the long range order is not destroyed. It was 
shown in Jl3| that the winding modes couple to the moments of the SU(N) magnetic flux 



associated with the gauge field A. 
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2.2. Chiral quantization 

The analysis of the matrix model simplifies considerably if it is formulated in terms 
of the chiral variables |5| 

- v; 

representing N x N hermitian matrices. The Hamiltonian in the new variables is 

H = ~Tr(X + X_ + X_X+), (2.8) 
where the matrix operators X± obey the canonical commutation relation 

[(X+Tj, = -i 5\5l (2.9) 

The partition function for fixed iV is now given by a path integral with respect to the one- 
dimensional hermitian matrix fields X + (x), X_(x) and A(x), satisfying periodic boundary 
conditions. The action (12.51) reads, in terms of the new fields, 



S = Tr f (iX+V A X- - \x + X_ ) dx. (2.10) 



o 



The partition function ( |2.10| ) depends on the gauge field only through the global holonomy 
factor, given by the unitary matrix 

VL = feC RA{x)dx . (2.11) 

This can be seen by imposing the gauge A = 0, which can be done at the expense of 
replacing the periodic boundary condition X±(f3) = X±(0) by a 5 , L r (A'")-twisted one 

X + (f3) = X + (0), X_((3) = tt^X-iO) Q. (2.12) 

2.3. Gauge-invariant collective excitations: momentum and winding modes 

The momentum modes V n jR of the matrix field are gauge-invariant operators with 
time evolution V n /n(t) = e~ nt / R V n /R(0). Here n can be any positive or negative inte- 
ger. The spectrum of momenta is determined by the condition that the time dependence 
is periodic in the imaginary direction with period 2nR. The operators satisfying these 
conditions are 

T/ _ p ±inx/R m y (^\n/R _ J TrX™ R if 71 > , . 

V n/R - e rrX ± {x) -< ]n]/R • (2.13) 



This form of the left and right momentum modes was first suggested to our knowledge by 
Jevicki 0. 

As it was shown in ||13|| , the winding modes V n R are associated with the SU(N) 



magnetic flux through the compactified spacetime 

V nR = Tr fl n (2.14) 
where O is the U(N) holonomy ( 2.1 1| ) . From the world sheet point of view the winding op- 



erator V n R creates a puncture with a Kosterlitz-Thouless vortex, i.e. a line of discontinuity 
2%Rn starting at the puncture. 
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2-4- The partition function as a three-matrix integral 



In the A = gauge the path integral with respect to the fields X±(x) is gaussian 
with the determinant of the quadratic form equal to one. Therefore it is reduced to the 
integral with respect to the initial values X± = X±(Q) of the action ( |2.10|) evaluated along 
the classical trajectories, which satisfy the twisted periodic boundary condition ( |2.12| ). 
Therefore the canonical partition function of the matrix model can be reformulated as an 
ordinary matrix integral with respect to the two hermitian matrices X+ and X- , and the 
unitary matrix Q: 



Z 



N 



J DX + DX_dQ e ^(x + x_- q x_nx + n-i)^ (215) 



where we denoted 

q = e ip = e 2lTlR . (2.16) 

A general perturbations by momentum modes can be introduced by changing the 
homogeneous measures DX± to 

[DX±] = DX± e ±lTrU±{x± \ (2.17) 

where we introduced the matrix potentials 

U±(X±) = Rj2t± n TrX!t /R (2.18) 

n>0 

Similarly, a general perturbation by winding modes can be introduced by changing the 
invariant measure DQ on the group U(N) to 

[DO] = Dfl e u( q 1/2 ^-u( q -^n) (2 1Q) 

where we introduced the following matrix potential 

U(Q) = ^t n Tr IT. (2.20) 

In case of a generic perturbation with both momentum and winding modes, the three- 
matrix model ( |2.15|) is not integrable. However, it can be solved exactly in case of an 
arbitrary perturbation only by momentum or only by winding modes. In the following we 
will discuss in details these two integrable cases. In both cases it is possible to integrate 
with respect to the angles and the matrix partition function reduces to an eigenvalue 
integral. 
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2.5. Integration with respect to the angles 



Let us consider the case of a perturbation with only momentum modes. Then the 
gauge field plays the role of a Lagrange multiplier for the condition = and the 

two matrices can be simultaneously diagonalized. Applying twice the Harish-Chandra- 
Itzykson-Zuber formula we can reduce the matrix integral to an integral over the eigenval- 
of the two hermitian matrices matrices X±. We write, after rescaling the 
integration variables, 

°? n 

Z N {t)= / \{[dxt][dx- k ] det(>j4) det(e-^ 7 <) (2.21) 

h— 1 

-OO 1 

where the measures are defined by 

[dx±] = dx± e zU±{x±) . (2.22) 
Then the grand canonical partition function can be written as a Fredholm determinant 

Z(fi, t) = Det(l + e-^K+KJ) (2.23) 

where 

[K+f\(x-)= [[dx+] e*»+*-f(x + ) 

J (2.24) 
[K_f](x+)= {dx_} e-**+*-f(xJ). 



3. The partition function of the non-perturbed theory 

3.1. Eigenf unctions of the non-perturbed Hamiltonian. 

In the absence of perturbation, the partition function is given by the Fredholm deter- 
minant 

Z(fj) = Det(l + e-^ +Ho) ) (3.1) 

and can be interpreted as the grand canonical finite-temperature partition function for 
a system of non-interacting fermions in the inverse gaussian potential. The Fredholm 
determinant ( |3.1|) can be computed once we know a complete set of eigenfunctions for the 
one-particle Hamiltonian 

H = --(x+x- + x-x + ). (3.2) 
The canonically conjugated operators x + and x_ 

— X-X + = —i (3.3) 
are represented in the space of functions of x+ as 

X— i&x+i 
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and in the space of functions of x — as 



X- = X-, x + = —id : 



Since the inverse gaussian potential is bottomless, the spectrum of the Hamiltonian is 
continuous 

H iP*(x±) = Ti(x±d x± + l/2)^f(x±) = £tyf (x±) (E e R). (3.4) 
The eigenfunctions in the x± -representation, which we denote by Dirac brackets 

^|(x±) = (E\x±) 

are given by 

(E\x±) = ^= x f E -^ (3.5) 

'TV 



where the phase factor 0o = (po(E) will be determined in a momenta. The solutions ( |3.5| ) 
have a branch point at x± = and are defined unambiguously only on the positive real 
axis. For each energy E there are two solutions, since there are two ways to define the 
analytic continuation to the negative axis. The wave functions relevant for our problem 
have large negative energies and are supported, up to exponentially small terms due to 
tunnelling phenomena, by the positive axisi. Therefore we define the analytic continuation 
to negative axis as 



(3.6) 



(E\ -x+) = (E\e-™x+) = e KE (E\x+) (x+ > 0) 
(E\ = (E\e i7r x_) = e 7zE (E\x-) (x_ > 0). 

The functions ( |3.5| ) with E real form two orthonormal sets 

oo 

J dx±(E\x±)(x±\E') =5(E -E'). 

— oo 

We will impose also the bi-orthogonality condition 

oo 

dx+dx- elx+x _ ^ = _ E , (3 g) 



(3.7) 



. 2tt 

oo 



2 The eigenfunctions ( |3.5| ) play in the chiral quantization the same role as the parabolic cylinder 
functions in the standard quantization, based on the original Hamiltonian (|2.1|). 



This corresponds to the theory of type I of ||. 
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which fixes the phase 0o 



e **o(i3) = yX e -f(s-i/2) r ( iS + 1/2 ). (3.9) 



The phase defined by ( |3.9| ) has in fact a small imaginary part, which we neglecto. This 
is because the theory restricted to the Hilbert space spanned on these states is strictly 
speaking not unitary due to the tunneling phenomena across the top of the potential. In 
the following we will systematically neglect the exponentially small terms 0(e ). With 
this accuracy we can write the completeness condition 

oo 

/ dE(x+\E)(E\x-) = --L e - ix + x -, (3.10) 
J v 2tv 



where the right-hand side is the kernel of the inverse Fourier transformation relating the 
x + and X- representations. For real energy the phase 4>o is real (up to non-perturbative 
terms), but we will analytically continue it to the whole complex plane. In this case the 
reality condition means 

ME) = ME). (3.11) 

3.2. Cut-off prescription, density of states and free energy 

To find the density of states, we need to introduce a cutoff A such that A >> \i. 
This can be done by putting a completely reflecting wall at distance = a/2A from 

the origin. Since the wall is completely reflecting, there is no flow of momentum through 
it, x+ ^~ = 0. Thus the cutoff wall is equivalent to the following boundary condition at 

x+ = x- = vA 

V^(VA) = V_(v / A) (3.12) 

on the wavefunctions Q3.5| ) . This condition is satisfied for a discrete set of energies E n ( n G 
Z) defined by 

0o(-B n )-S n logA + 27rn = O. (3.13) 
From ( |3.13j ) we can find the density of the energy levels in the confined system 

Now we can calculate free energy !F(fi,R) = logZ(fi, R) as 

/oo 
dEp(E) log 1 + e-^+E) ( 315 ) 
- oo 



The pure phase factor would be ^J cos ^ E Y{iE + 1/2) = yT + e* E e 



i<t> (E) 
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with the density ( |3.14| ). Integrating by parts in (|3.15| ) and dropping out the A-dependent 
piece, we get 

m = ~J dME) log (l + e-^+*)) = -R d E i ^ +E) . (3.16) 

We close the contour of integration in the upper half plane and take the integral as a sum 
of residues. This gives for the free energy 

T = -i <t>o(ir/R-n)- (3.17) 

r=n+i>0 

Using the explicit form ( ^.9| ) of <fio we can represent the free energy as a sum over pairs of 
positive half-integers 

oo 

T{^R) = log (/j, + ir + is/ R) (3.18) 

r,s>l/2 

which is explicitly invariant under the T-duality R — > 1/R, \i — > Rfi. From here it follows 
that the free energy satisfies the functional equation 

4:Sm{dj2R)sm(d t ,/2)J 7 (fi,R) =-log/x. (3.19) 

From (|3.17| ) we can express the phase factor in terms of the free energy 

<£o = -2sin (d fl /2R)J r (fi 1 R). (3.20) 

3.3. Calculation of the free energy using the E -representation of the operators x± 

Let us now give an alternative, algebraic derivation of the free energy of the non- 
perturbed theory. For this purpose we will write the canonical commutation relation ( |3.3j ) 
for the matrix elements of the operators x± in ^-representation. 

The action of the operators x± to the wave functions is equivalent to a shift of the 
energy by an imaginary unit (we consider the wave functions as analytic functions of E) 
and a multiplication by a phase factor: 

(E\x±\x±) = e ±HMETi)-Mm ( E T i\ x± } (3.21) 

As a consequence, the operators x± are represented by the finite difference operators acting 
in the -E-space 

x± -> e^e^ dE e ± ~^. (3.22) 
The Heisenberg relation (|3.3| ) is equivalent to the condition 

ME + \i) - ME -\i) = - log(-S) (3.23) 



which is equivalent to the functional constraint (|3.19| ). 
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4. Integrable perturbations by momentum modes. 

4-1. One-particle eigenf unctions and the density of states of deformed system 

Now we will consider perturbations generated by momentum modes V n m- As we 
argued in sect. 2.4, this can be achieved by deforming the integration measures dx± to 

[dx±] = dx± exp (±iU±(x±)) , U±(x±) = R^^t±kX± . (4-1) 

fc>i 

In this section we will show that such a deformation is exactly solvable, being generated 
by a system of commuting flows H n associated with the coupling constants t± n . The 
associated integrable structure is that of a constrained Toda Lattice hierarchy. The method 
is very similar to the standard Lax formalism of Toda theory, but we will not assume that 
the reader is familiar with this subject. It is based on the possibility to describe the 
perturbation by vertex operators as deformations of the ^-representation of the canonical 
commutation relation. 

The partition function of the perturbed system is given by the Fredholm determinant 
( j2.23|) , where the integration kernels ( |2.23| ) and ( |2.24j ) are defined with the deformed mea- 
sures Q4.1|) . The deformed kernel is diagonalized by a complete orthonormal system of 
wave functions! ^±(x±) 

J dx±V%(x±)V%(x±) = 5(E - E') (4.2) 
which satisfy the following three conditions: 

1) they are eigenfunctions of the evolution operator relating the points x = and x = (3 

e -PH ^E = e -PE^ (43) 



2) they behave at infinity as 

#f (x±) ~ xf E ~ h e * (E ) e ^±(*±), (4.4) 



3) they satisfy the bi-orthogonality condition 



dx+dx- qE( x _) e i3: + x -^'(x + ) = 6(E-E / ). (4.5) 



2,71 



The new eigenfunctions are not necessarily eigenvalues of the Hamiltonian Hq itself. 
Indeed, the evolution operator e~^ H ° acts trivially on any entire function of x± R '. In 
particular, the condition (|4.4|) is compatible with the condition (|4.3|) . 



5 As before, we understand the completeness in a weak sense, i.e. up to non-perturbative terms 
due to tunnelings phenomena, which we have neglected. 
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Once the basis of the perturbed wave functions is found, the logarithm of the Fredholm 
determinant can be calculated as the integral ( |3.15|) , where the new density of states is 
obtained from the phase 4>(E) in the asymptotics (fOj) 



Therefore the phase (f>(E) contains all the information about the perturbed system. It is 
related to the free energy by 

T= -i <Kir/R-fj) ■ (4.7) 

r>l/2 

<f>{-y) = 2sin (hd^/2R)F(^ R). (4.8) 

4-2. Dressing operators 

Now we proceed to the actual calculation of the phase (p. The method is a generaliza- 
tion of the algebraic method we have used in sect. 3.3. First we remark that the perturbed 
wave functions are given by the rhs of (|4.4|) up to factors of the form 



W±(x±) =exp ^iRj2 v ±n x± n/R j . (4.9) 
which satisfy the second condition (|4.3|) and tend to 1 when x± — > 00. The unknown 



coefficients v n ) and the phase <fi are determined functions of E and the couplings t n by the 
third condition (|4.5j ). 

It follows from (|3.22| ) that the deformed wave function 

ttf (x±) = e tU±{x±) e^ {E) x% %E ~ l2 W ± {x ± ) (4.10) 

can be obtained from the bare wave functions ( |3.5|) by acting with two finite-difference 
operators in the -E-space VV+ and VV_ 

*f(x±) = (E|e ± 5^°VV ± |x±). (4.11) 

The explicit form of the dressing operators VV± is obtained by replacing x± — > e^ ldE in 
( |4.10|) . The dressing operators are unitary 

= W-Wl = 1 (4.12) 

since they relate two orthonormal systems of functions. Further, the bi-orthogonality 
condition (|4.5| ) is equivalent to the identity 

wie^°W+ = 1. (4.13) 

The identity ( pp means that the product VV-VV.). 1 does not depend on the perturbation, 
which is a general property of all Toda lattice systems ||. 
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4-3. Lax operators and string equation 



We have seen that the partition function is expressed in terms of the phase 0(— //). 
Therefore we assume that the energy is at the Fermi level E = — //, and consider all 
variables as functions of \x instead of E. Let us denote by cD the operator 



Gi = e id » (4.14) 



shifting the variable fx by i. The operators cD and \i satisfy the Heisenberg-Weyl commu- 
tation relation 

[u,fi]=u, [u>~ , fj] = — Cj~ . (4-15) 

Now let us consider the representation of these commutation relations in the perturbed 
theory. The dressing operators W± are now exponents of series in cD with ^/-dependent 
coefficients 

Vv ± =e ^^»>i t± """ /fl e =F^(M) e^E^^i-M -"" /fl . (4 . 16) 

The operators 

L+ = W+wW7 1 , L_=W_w- 1 Wl 1 , 

known as Lax and Orlov-Schulman operators satisfy the same commutation relations as 
the operators Cj and \x 

[L+, M+] = iL+, [L-, M_] = -%L— (4.18) 

The Lax operators L± represent the canonical coordinates x± in the basis of perturbed 
wave functions, 

{E\e ± ^ i4)a W ± L ± \x ± ) = {E\e ± ^ i<t,0 W ± x ± \x ± ) (4.19) 

while the Orlov-Shulman operators M± represent Hamiltonian Hq = — |(x_|_x_ + £_£+). 
Therefore the L and M operators are related also by 

M + = M- = -(L + L- + L-L + ). (4.20) 

The last identity is not satisfied automatically in the Toda system and represent an addi- 
tional constraint analogous to the string equations in the minimal models of 2D quantum 



gravity. The relation ( 4.20 ) is proved by inserting the identity ( 4.13 ) in the bare relation 
for E = —fi, satisfied by the operators ( |3.22|) . The string equation can be also written as 
the Heisenberg commutation relation between the two Lax operators 

[L+,L_] = -i (4.21) 

which also follows directly from (|3.3|). 
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The operators M± can be expanded as infinite series of the L-operators. Indeed, as 
they act to the dressed wave functions as 



(E\ e ± ^° W± M±\x±) = ±i(x±d x± - l/2)*|(x±) 

= ( kt ±* x ± R + » + E v ±* x ± k/R l *± ( x ±)- 



(4.22) 



,fc>i fc>i 
we can write 

k/R 



M± = J2 kt ±kL± R + V + v ±kL± /R . (4.23) 



fc>i fe>i 



In order to exploit the Lax equations (|4 . 1 8| ) and the string equations (|4.20| ) we need the 
explicit form of the two Lax operators. It follows from fl4.16| ) that L± can be represented 
as series of the form 



v + = e ^' \uj + 2_^ a k UJ \ e 

k>i 



(4.24) 



fc>i 



Integrable flows 

Let us identify the integrable flows associated with the coupling constants t n . From 
the definition (|4.17[ ) we have 

d tn L± = [H nj L±], (4.25) 
where the operators H n are related to the dressing operators as 

H n = (d tn W+)W+ 1 = (a t „W_)Wl 1 . (4.26) 

The two representations of the flows H n are equivalent by virtue of the relations ( |4.12|) and 
( |4.13| ) . A more explicit expression in terms of the Lax operators is derived by the following 
standard argument. Let us consider the case n > 0. From the explicit form of the dressing 
operators it is clear that H n = W+u n / R W+ 1 + negative powers of uj l t R . The variation of 
t n will change only the coefficients of the expansions ( |4.24| ) of the Lax operators, preserving 
their general form. But it is clear that if the expansion of H n contained negative powers 
of cD 1 /^, its commutator with L_ would create extra powers u~ 1 ~ k t R . Therefore 

H ±n = {L n J R )> + \{L n J R ) Q , n>0, (4.27) 

where the symbol ( ) > means the positive (negative) parts of the series in the shift operator 

Cj 1 / r and ( )q means the constant part. By a similar argument one shows that the Lax 
equations ( |4.25[ ) are equivalent to the zero-curvature conditions 

d tm H n - d tn H m - [H m , H n ] = 0. (4.28) 
13 



Equations ( fl25|) and ( |4.25|) imply that the perturbed theory possesses the Toda lattice 
integrable structure. The Toda structure implies an infinite hierarchy of PDE's for the 
coefficients v n of the dressing operators, the first of which is the Toda equation for the 
phase 4>(n) = <fi(E = —fi) 

i - - <p(fi) = e i ^t J -)- i ^t J '- i / R ) — e i4>(p+i/R)-i4>(p) ^ (4.29) 

The uniqueness of the solution is assured by appropriate boundary conditions |T3|], which 
are equivalent to the constraint ( f4.21| ). 

4-5. Representation in terms of a bosonic field 

The momentum modes can be described as the oscillator modes of a bosonic field 
<p(x+, X-) = (p+(x+) + (f-(x-). The bosonization formula is 

^l = ~^~ l/2 {x ± ) = Z~ x e ±lv±{x±) ■ Z (4.30) 
where Z is the partition function and 

<p±(x±) = +R^2t k x^/ R + + ^iogx ± - T x ± k/R 7fr- ( 4 - 31 ^ 

k>i k>i k 

Then by Q4.22]) the operators M± are represented by the currents x±d±(f 

Mt^(x ± )\ E= _^ l/2 = Z- 1 x ± d ± ip ■ Z. (4.32) 

4-6. The dispersionless (quasiclassical) limit 

Let us reintroduce the Planck constant by replacing \i — > n/% and consider the quasi- 
classical limit h — ► 0. In this limit the integrable structure described above reduces to the 
dispersionless Toda hierarchy |1^,|0|,^[ , where the operators \x and u can be considered as 
a pair of classical canonical variables with Poisson bracket 

{uj,fi} = uj. (4.33) 

Similarly, all operators become c-functions of these variables. The Lax operators can be 
identified with the classical phase space coordinates x±, which satisfy 

{x+,x_} = l. (4.34) 

The two functions x±(u, E) define the classical phase-space trajectories as functions of the 
proper time variable r = \oguj. 

The shape of the Fermi sea is determined by the classical trajectory corresponding to 
the Fermi level E = —fx. In the non-perturbed system the classical trajectory is 

x+(u) = y/jl lo, x-(u) = yfjl uj~ 1 (4.35) 
14 



and the Fermi sea has a hyperbolic shape 

x+x- = /i. (4.36) 
In the perturbed theory the classical trajectories are of the form 

x ± =L ± {u,fj,) (4.37) 
where the functions L± are of the form 

L±{u,n) = e* 9 »+u ±1 ^l + J2a ± Mu^ k / R ^j . (4.38) 

The flows H n become Hamiltonians for the evolution along the 'times' t n . The uni- 
tary operators W± becomes a pair of canonical transformations between the variables 
u>, \i and L±,M±. Their generating functions are given by the expectation values 
S± = Z~ x (fi±(x±) ■ Z of the chiral components of the bosonic field <p 

S± = ±RJ2 t±k x k J R +^logx ± -^±Rj2ln xt IR (4.39) 
k>i k>i 

where Vk = dT jdtk- The differential of the function S± is 

dS± = M±dlogx± + logiv dfi + R^H n dt n . (4.40) 

If we consider the coordinate w as a function of either x+ or x_, then 
The classical string equation 

x + x- = M + = M- (4.42) 
can be written, using the expansion (|4.23|) of M±, as 



x + x- = kt k x+ R + fx + ^2 v k ^+ 

k>l k>l 

x+X- = kt-k x k J R + fx + V-k x_ k 
k>i k>i 



(4.43) 



The first of these expansions is convergent for sufficiently large and the second one for 
sufficiently large x_. Comparing the two equations one can extract the form of the Fermi 
surface. Technically this is done as follows |17]. First, note that if all t±k with k > n 



vanish, the sum in ( |4.38| ) can be restricted to k < n. Then it is enough to substitute the 
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expressions ( [4.38 ) in the profile equations ( 4.43 ) and compare the coefficients in front of 
The two expansions ( f4.43| ) can be combined into a single equation 

x-x- - ^2 kt k H k (u) = \l (4.44) 

fc#0 

where 

H ±k (u) = [L k J R (u;)]> + ±[L k J R [u>)] (k > 0). (4.45) 

The left hand side can be interpreted as the Hamiltonian for the perturbed system. It 
defines the profile of the perturbed Fermi sea, which is a deformation of the hyperbole 

(ill). 

4-7. One-point correlators in the dispersionless limit 

It follows from logw = d fJi S±(x±) that the ^/-derivative of the one-point correlators 

{TrX n /R) = dT (/U6) 

is given by the contour integrals 

d 2 T 



d\i dt± 7 



— <tdu l,R [L±(u)} lnl/R (n > 1) (4.47) 

TCI J 



where the closed contour of integration in the variable uj x I r goes along the arc between 
e -nrR an( j w _ gJTri? (N t e that the integrand is expanded in Laurent series in uj 1 ^.) 



to 



4-8. Example: sine- Gordon field coupled to 2D gravity 

The simplest nontrivial string theory with time-dependent background is the sine- 
Gordon theory coupled to gravity known also as Sine-Liouville theory. It is obtained by 
perturbing with the lowest couplings t\ and t—\. In this case 



x± 



e^V^l + aiiw^). (4.48) 



and (fOS) give 



\i e 



a ±1 =t Tl e-^-^t 
The first equation is an algebraic equation for the the susceptibility 

u = d 2 ^T = -Rd^. 



(4.49) 



It was first found (for the T-dual theory) in |13|]. This algebraic equation resumes the 
perturbative expansion found in ||21|| . The equation ( 4.48 ) with a± given in Q4.49|) was first 
found in |16|] by integrating the Hirota equations for the Toda hierarchy with the boundary 
condition given by the non-perturbed free energy ( |3.18| ). 
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5. Integrable perturbations by winding modes 

5.1. The partition function as integral over the U(N) group 

Let us consider a perturbation of the c = 1 string theory by only winding modes. 
The relevant variable in this case is the gauge field A(x) and the matrix model ( |2.4[ ) / can 
be viewed as a two-dimensional gauge theory defined on the disc having as a boundary 
the compactified Euclidean time interval. The only nontrivial degree of freedom in such a 



theory is the holonomy factor around the circle ( 2.11 ). The canonical partition function is 
given by the U(N) -integral 



Z 



N 



[DO] 



|g-V2 ft® J-gl/2 I X 0| 



Z N (t r 



A/2 



-1/2 1 -JV 



N 

n 

fc=i 



[<fefc] 



n 



_1 / 2 Zj — q 1 / 2 Zk 



where the integration goes along the unit circle \z\ = 1 with measure 



(5.1) 



where I is the unit N xN matrix. The integrand depends only on the eigenvalues Z\ , ... , Zn 
of the unitary matrix Q and the partition function becomes 



(5.2) 



[dz]=dze°^ z ^- G ^ x \ (5.3) 
This representation of the partition function was studied (for the non-perturbed theory) 



by Bulatov and Kazakov [12]. Note that the absolute value can be abandoned, since the 
integrand is homogeneous. 

The grand canonical partition function can be written as the absolute value of a 
Fredholm determinant 



Z(n,t r 



Det(l + q^K) 



where the Fredholm kernel is defined by the contour integral 

[dz] f(z') 



(Kf)(z) 



2%i q 1 / 2 z — q 1 / 2 z' 



(5.4) 



(5.5) 



The integral has poles when Zi = qzj and should be evaluated by adopting a prescription 
for surrounding the poles. The prescription used in [|T2[ is to add a small imaginary part 
to R so that \q\ < 1. In the case of real q e [0, 1) the partition function was studied by M. 



Gaudin [22 



5.2. Evaluation of the partition function of the non-perturbed theory 



Boulatov and Kazakov showed in |12| that the prescription for the contour integration 
gives for the Fredholm determinant (|5.4| ) the same result as (|3.18| ) , up to non-perturbative 
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terms 0(e 7rM ). Here we recall their calculation. The integration kernel (5J3) acts to the 
monomials z n as 

Kz n = {q n+ h" ifn>0 (56) 
lO ifn<0 
and the Fredholm determinant (|5.4|) reads 

r>0 

The free energy then is equal to an infinite sum 

/ \m n imfx i 

n>0 m>l y y 

The last sum is the result can be written as the integral 

= f ~ 7 ■ h( q Z y ■ u y (5-8) 
J y 4smh(yR) smh(y) 

c 

along a contour C circling around the poles y = in, n > only. If |g| < 1, then C goes 
from — oo to and then up the imaginary axis to zoo. 

Note that if we close the contour around the poles of smhirRy, then the result will be 
the sum over the residues k = inR (n > 0) which can be written as the partition function 
for the dual radius R = 1/R. 

5.3. Toda integrable structure 

The Fredholm determinant ( |5.4j) with non-homogeneous measure ( |5.3|) can be repre- 
sented as the Fock expectation value in a theory of chiral fermions defined on the unit 



circle |T3[]. The deformations by winding modes are introduced as Bogolyubov transfor- 
mations of the left and right fermionic vacua and the partition function was identified as 
a r-function of the Toda Lattice hierarchy. 

As a consequence, the free energy satisfies an infinite hierarchy of PDE with the Toda 
'times' t n . The first one is the Toda equation dual to the equation ( |L29|) 

— S—d>(u) = e ik»)-ik»-i) _ e #0+i)-W (5 9) 

dhdt-i J v ; 

where 

0(/x) = 2sin(a M /2)^(/x, J R). (5.10) 

Using the scaling relation given in the Appendix, one can reduce ( |5.9| ) to an ordinary 
differential equation, which should be solved with initial condition ( |5.8|) . In the quasi- 
classical (genus zero) limit, this differential equation can be integrated to an algebraic 
equation for the string susceptibility uq = d^cf) 

fie ^ d ^+ iJ-^R - l)e^ 9 ^ = 1. (5.11) 

The one- and two-point correlators were calculated by Alexandrov and Kazakov from the 
higher equations of the Toda hierarchy [|16[. These results were later confirmed in [|17 



using the Lax formalizm of the constrained Toda system. 
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5.4- Lax operators and string equation 

We will consider the potential U(z) in the measure ( |5.3[ ) as the value on the unit circle 
zz = 1 of the potential 

U(z, z) = + l -^ n ) ( 5 - 12 ) 

n>0 

defined in the whole complex plane. Let us assume that the spectral variables z and z are 
represented (in some sense) by a pair of Lax operators of the form 

V , " ' (5-13) 

L = e *** Co' 1 j 1 + u k Co l ~ k e - *'*, 

where the phase <p is related to the free energy by 

4>{n) = -i[.F(/x + i/2) - ^(A* - */2)]- (5-14) 

In order to find the constraint satisfied by L and L, we consider the non-perturbed 
theory, for which the expressions in the parentheses are equal to one. The functional 
equation ( |3.19| ) means that the bare phase (f>o satisfies 

<p (fj, + i/R) - o (a* _ = ~i l°g A*- (5.15) 
In the case of no perturbation, this is equivalent to the algebraic relations 

L l/R L 1/R + L l/R L 1/R = fi (516) 

and 

[L 1 ' R ,L 1 l R }=i. (5.17) 

The second identity is invariant with respect to the dressing procedure and therefore are 
satisfied by the the general Lax operators (|5.13 ). This is the string equation for the 
constrained Toda system. 

5.5. Quasiclassical limit and relation to the conformal map problem 

In the dispersionless limit — > oo the two Lax operators define a smooth closed curve 
7 in the complex plane, whose equation is written in a parametric form as 

z = L(u), z = L(u) (M = l). (5.18) 

For example, in the case when only t\ and t-i are nonzero, the explicit form of the curve 
7 is 

2 R . (5.19) 
z = e-* 9 ^ uj- 1 (1 + t ie — d ^ u) R . 
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Assuming that the couplings t n and t_ n are complex conjugate, the map ( |5.18|) can 
be extended to a conformal map from the exterior of the unit disk uj < 1 to a connected 
domain T> with the topology of a disk containing the point z = oo and bounded by the 
curve 7. The couplings t n (n 7^ 0) and fx can be thought of as a set coordinates in the 
space of closed curves. The relation between the conformal maps and the dispersionless 
Toda hierarchy was studied recently in a series of papers 



6. Conclusion 

In these notes we explained how the the perturbations of the compactified c = 1 
string theory by momentum or winding modes are described by integrable deformations of 
a gauged matrix model on a circle. The momentum and winding modes are associated with 
the collective excitations of the matter and gauge fields. The deformed system is described 
by a constrained Toda Lattice hierarchy. The Lax formalism for the Toda lattice hierarchy 
allows to calculate explicitly the free energy and the correlation functions of the electric 
or magnetic operators for any genus. In particular, the partition function is a r-function 
of the Toda lattice hierarchy. 

The integrability takes place only in the grand canonical ensemble, in which the size 
N is a dynamical variable and the string interaction constant is controlled by the chemical 
potential \x. The partition function is a Fredholm determinant, and not a usual determinant 
as in the case of the open matrix chains for the c < 1 string theories. 

It is not likely that the integrability is conserved for perturbations with both moment 
and winding modes. Nevertheless, the calculation of the correlation functions of the mo- 
mentum modes in presence of a perturbation modes (or vice- versa) seem to be performable 
albeit very difficult. This calculation might help to check the hypothesis (related to the 
FZZ conjecture [[Kj) that a strong perturbation by winding modes can lead to a curved 
background with Euclidean horizon |13| . 



Another problem, also related to the conjecture, is to find the limit \ijt\t-\ — > 



of the matrix model, which is relevant to the sine-Liouville string theory. This limit seems 
to be subtle because at small \i the non-perturbative effects can enter into the game. In 
any case, the \x — > limit of the correlators of the matrix model does not seem to reproduce 
the results obtained in the Sine-Liouville theory []T0|j2l 



Finally, it would be very interesting to understand the origin of the integrability from 
the point of view of the world-sheet string theory. 
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Appendix A. Momentum and winding modes in the c = 1 Euclidean string 
theory 

In this section we will recall briefly the world-sheet description of the c = 1 Euclidean 
string theory with compact target space. The elementary excitations in this theory are 
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closed surfaces, or string world sheets, embedded in a circle of radius R. An embedded 
surface is defined by metric tensor g a b(cr) (a, b = 1, 2) and the position x(a) in the target 
space as functions of the local coordinates a = {01,(72}. The free energy of the string 
theory is given by the functional integral over all connected surfaces 

^F(R, fi,h) = J Vg ab Vx e~ s ^) (A.l) 

with a weight given by the the Polyakov action 

S(g ab ,x) = -^ J d 2 a y /d^tg'[g ab d a xd b x + 4:7rii + TZW\ogh}, (A.2) 

world sheet 

where TZ is the local curvature associated with the metric g ab . The parameter \i is the 
cosmological constant, coupled to the area of the world sheet, and K is the string coupling 
constant, which is associated with the processes of splitting and joining of closed strings. 
By the Euler theorem the global curvature is related to the genus h of the world sheet as 

-!- J d 2 (J^/fetg~ = 2 - 2h (A.3) 

and one can write the free energy as a series 

F(R,g) = J2h 2h F {h) (R^)- (A.4) 

h>0 

In the conformal gauge g ab = e~ 2 ^^8 ab , the conformal factor (j> becomes a dynamical 
field due to the conformal anomaly, and the world-sheet action becomes essentially a c = 1 
conformal field theory coupled to a Liouville field 

S=^- [ d 2 a [(dxf + (d(j)) 2 + /je b ^ + n {2) (logh + Q(j))+ ghosts]. (A.5) 

47T J 

world sheet 

The background charge Q and the exponent b of the Liouville field are determined by the 
requirement that the total conformal anomaly vanishes 

Cmatter + CLiouville + Cghosts = 1 + (1 + QQ 2 ) - 26 = (A. 6) 

and that the perturbation due to the cosmological term fxe b( ^ is marginal 

6 ^- 6 ' = i. (A.7) 
4 

These two conditions give 

Q = -2, b=-2. (A.8) 
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(The value Q = 2 does not lead to a sensible classical limit.) The invariance of the action 
with respect to shifts <fi — > (p + </>o implies that the free energy depends on \i and % through 
the dimensionless combination \ih 1 which is the statement of the double scaling limit in the 
c = 1 string theory. Therefore we are free to choose % = 1 and write the genus expansion 
QA.4 ) as an expansion in 1 / /j, 2 . 

The primary operators associated with the matter field x(a) are the vertex operators 
V e (a) ~ e -iex(a) an( j £h e Kosterlitz-Thouless vortices V m (a). The operator V m (a) is 
associated with a discontinuity lixm of the field x around the point a on the world sheet. We 
call e and m electric and magnetic charges, in analogy with the Coulomb gas on the plane. 
The electric and magnetic charges should satisfy the Dirac condition em = 2ttx integer. 
From the point of view of the compactified string theory observables, the electric and 
magnetic charges are the momentum and winding numbers, correspondingly. The spectrum 
of the electric and magnetic charges for a periodic target space x + 2nR = x is 

e = n/R, m = nR (n £ Z). (A. 9) 

If we write the position field i as a sum of a holomorphic and antiholomorphic parts, 
x = xr+xl, then the Kosterlitz-Thouless vortices are described by the vertex operators for 
the dual field x = xr — xl, whose target space is the circle of radius 1/R. When integrated 
over the world sheet, these operators should be accompanied by nontrivial Liouville factors 
compensating their anomalous dimensions. The integrated operators have the form 

V e ~ f rfVe-"e (e - 2) ^, 

J (A. 10) 

The Liouville exponents are determined by the condition that the integrands are densities. 

We are interested in deformations of the string theory obtained by allowing electric 
charges e = n/R with fugacities t n and magnetic charges m = nR with fugacities t n . This 
is achieved by adding to the action ( |A.5|) the perturbation term 

5S = ^2(t n V n/R + t n V nR ). (A. 11) 

The translational invariance of the functional measure T>cj) yields the following Word 
identity for the free energy 



~ 2 "f + E M " ^) + E {I - 2) f«f - 0. (A,2) 



This means that the couplings t n , t n and the string coupling h scale with respect to the 
cosmological coupling [i 

tn~/^ H/i? , tn-M 1 -^ 1 , (A.13) 
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The T-duality symmetry of the original theory 



x^x, R^l/R, n^n/R (A.14) 

holds for the perturbed theory if one also exchanges the couplings as t n <-> t n (up to a 
rescaling by an R- dependent factor). 

In the perturbation (|A.11|) we should retain only the relevant charges |e| < 2/R and 
\m\ < 2R. However, the correlation functions of a finite number of irrelevant operators are 
perfectly meaningful. The interesting phase of the deformed theory, which can be thought 
of as Coulomb gas coupled to quantum gravity, is the plasma phase, where the fugacities 
of the charges are tuned so that the Debye length is of the order of the size of the 2D 
universe. 
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